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THERMODYNAMIC PROPERTIES OF POTASSIUM TO 2100° K 

by Sheldon Heimel 
Lewis Research Center 

SUMMARY 

Consistent thermodynamic properties of potassium from the melting point 
(336. 35° K) to 2100° K were computed for both the ideal and the real species and are 
presented in tabular form. 

The properties computed for the real monomeric vapor on the saturation line are as 
follows: pressure, volume, heat capacity, enthalpy of vaporization and enthalpy, entropy 
of vaporization and entropy, fugacity- pressure ratio, virial coefficients (second, third, 
and fourth), compressibility factor, density, partial derivative of volume with respect to 
temperature at constant pressure, and partial derivative of volume with respect to pres- 
sure at constant temperature. 

The properties computed for the ideal species (condensed phases, monomeric gas, 
and dimeric gas) are heat capacity, enthalpy, entropy, and free energy. 

The real-gas properties were computed from a virial equation of state. The equation 
used existing second virial coefficients based on the Rydberg potential. Third and fourth 
virial coefficients exhibiting limited aspects of hard-sphere behavior were adjusted to ex- 
perimental pressure-volume-temperature (PVT) data. 


INTRODUCTION 

Three areas requiring a reliable equation of state of the alkali metals are 
(1) Rankine- cycle engines, (2) radiation heat transfer in plasmas, and (3) calculation of 
transport properties. In one study of the Rankine-cycle engine at the Lewis Research 
Center (ref. 1), an equation of state for potassium was used to calculate the critical 
weightflow rate in a nozzle. 

There are two recent virial equations of state for potassium (refs. 2 to 5). These 
equations yield discrepant sets of thermodynamic properties. Therefore, the present 
work was performed to eliminate the discrepancies by presenting an improved equation. 
The improvement lies in the representation of each virial coefficient (up to the fourth) by 
the ’’best" available potential model. 



Second virial coefficients derived from the Rydberg potential (1100° to 1600° K) were 
taken from reference 6 . The dependencies of the third and fourth virial coefficients on 
the second virial coefficient were taken from the hard- sphere model (ref. 7). These de- 
pendencies correlate the experimental volumes (1127° to 1575° K) (refs. 2 to 4) better 
than analogous dependencies from the Lennard-Jones 12-6 potential (ref. 8 ). 

Analytic functions approximating the virial coefficients were used in the range from 
336. 35° to 2100° K in order to extend the range of the equation of state beyond the range 
of the pressure- volume-temperature (PVT) data. In addition, smoothed vapor pressures 
over the entire experimental range of 372° to 2078° K (refs. 2 to 4, 9 to 11, and unpub- 
lished data of K. J. Bowles of Lewis) were used. The justification for this procedure is 
the reasonable trend shown by calculated heat capacities, fugacity-pressure ratios, and 
compressibility factors. 

The results of the analysis are as follows: 

(1) An improved equation of state is obtained for the temperature range 336. 35° 
to 2100° K. 

(2) Volumes, heats of vaporization, and isentropic expansion indexes derived from 
the equation are in good agreement with experimental data. 

Tables of thermodynamic functions are presented for the ideal crystal, liquid, mo- 
nomer, and dimer and for the real monomer on the saturation line. 

THERMODYNAMIC PROPERTIES OF IDEAL SPECIES 

The thermodynamic properties of the condensed phases (crystal and liquid) and of 
the monomeric and dimeric gases are given in tables I, n, and HI, respectively. These 
properties were obtained in the following manner. 


Crystal 

The JANAF data (ref. 12) for the heat capacity at constant pressure for the standard 
state C°, the sensible enthalpy for the standard state - Hq, the entropy for the stan- 
dard state S° from 100° to 300° K were extrapolated to the melting point by the 
Lagrangian formula. (Superscript letter o indicates the ideal state. All symbols are 
defined in appendix A.) The melting point (336. 35° K) and the heat of fusion of 558 cal- 
ories per mole (2334.67 J/mole) were taken from Douglas et al. (ref. 13). 

The arbitrary base for assigning values to the enthalpy of crystal and liquid potas- 
sium (K) was a value of zero at 298. 15° K (i. e. , (h^qq 15 ) = o). Since table I gives 

( H 298. 15 " H o) as 1693. 0 calories per mole (7083. 5 J/mole), (Hg) = -1693. 0 calories 
^ * u/ c 'c 

per mole (-7083. 5 J/mole). 
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TABLE I. - THERMODYNAMIC PROPERTIES OF CRYSTAL AND LIQUID POTASSIUM 


Temper- 

ature, 

T, 

°K 

Heat cap- 
acity for 
standard 
state at 
constant 
pressure, 

C p/ R 

Sensible 
enthalpy 
for stand- 
ard state 
at T° K, 

h t “ H o 

Entropy 

for 

standard 

state, 

S°/R 

Sensible 
free en- 
ergy for 
standard 
state at 
T° K, 

-( f t - H o) 

Sum of 
sensible 
enthalpy 
at T°K 
and 

chemical 
energy at 
0° K for 
standard 
state, 

h£/rt 

Formation from assigned 
reference state 

Enthalpy 
change 
at T°K, 

( ah t)/ rt 

Log of 
equilib- 
rium 
constant, 
log Sf 

RT 

RT 

100 

2.9625 

2. 0834 

4. 2820 

2. 1986 

-6.4363 

0 

0 

200 

3. 2478 

2.6017 

6.4338 

3.8321 

-1.6581 

0 

0 

298. 15 

3. 5478 

2.8575 

7.7784 

4. 9209 

.0000 

0 

0 

300 

3. 5568 

2.8617 

7.8005 

4. 9388 

.0218 

0 

0 

a 336. 35 

3. 9624 

2. 9594 

8.2303 

5. 2709 

.4264 

0 

0 

a 336. 35 

3.8655 

3.7942 

9.0651 

5. 2709 

1.2612 

0 

0 

400 

3.7885 

3.7993 

9.7284 

5. 9292 

1.6694 

0 

0 

500 

3.6918 

3.7870 

10. 5629 

6.7759 

2.0830 

0 

0 

600 

3.6246 

3.7651 

11.2296 

7.4645 

2.3452 

0 

0 

700 

3. 5870 

3.7420 

11.7851 

8.0431 

2. 5249 

0 

0 

800 

3. 5790 

3.7218 

12.2633 

8. 5414 

2.6569 

0 

0 

900 

3.6005 

3.7069 

12.6858 

8.9789 

2.7602 

0 

0 

1000 

3.6516 

3.6985 

13.0675 

9.3690 

2. 8466 

0 

0 

b 1030. 90 

3.6734 

3.6975 ' 

13. 1789 

9.4815 

2.8710 

0 

0 

1100 

3.7323 

3.6977 

13.4191 

9.7214 

2.9232 

0 

0 

1200 

3.8425 

3.7050 

13.7483 

10. 0434 

2.9950 

0 

0 

1300 

3. 9823 

3.7208 

14.0612 

10.3405 

3.0654 

0 

0 

1400 

4. 1517 

3.7453 

14. 3624 

10.6171 

3. 1368 

0 

0 

1500 

4.3507 

3.7789 

14.6554 

10.8765 

3.2109 

0 

0 

1600 

4. 5792 

3.8216 

14. 9433 

11. 1217 

3.2891 

0 

0 

1700 

4. 8372 

3. 8736 

15. 2286 

11.3549 

3.3725 

0 

0 

1800 

5. 1249 

3.9350 

15. 5130 

11. 57 80 

3.4617 

0 

0 

1900 

5. 4421 

4. 0058 

15.7985 

11.7927 

3. 5574 

0 

0 

2000 

5.7888 

4. 0862 

16.0863 

12.0001 

3.6602 

0 

0 

2100 

6. 16 52 

4. 1761 

16.3777 

12. 2016 

3.7704 

0 

0 


a Melting point. 
^Normal boiling point. 



TABLE n. - THERMODYNAMIC PROPERTIES OF IDEAL POTASSIUM GAS (Kj) 


Temper- 

ature, 

T, 

°K 

Heat cap- 
acity for 
standard 
state at 
constant 
pressure, 

c°/R 

Sensible 
enthalpy 
for stand- 
ard state 
at T°K, 

H T ■ H o 
RT 

Entropy 

for 

standard 

state, 

S°/R 

Sensible 
free en- 
ergy for 
standard 
state at 
T° K, 

-( f t ' H o) 

Sum of 
sensible 
enthalpy 
at T°K, 
and 

chemical 
energy at 
0° K for 
standard 
state, 
H° /RT 

Formation from assigned 
reference state 

Enthalpy 
change 
at T°K, 
(AH°) f /RT 

Log of 
equilib- 
rium 
constant, 
log /T f 

RT 

100 

2. 5000 

2. 5000 

16.5404 

14.0404 

101.6531 

108.0894 

-41.6189 

200 

2. 5000 

2. 5000 

18.2732 

15.7732 

52. 0766 

53.7347 

-18. 1949 

298. 15 

2. 5000 

2. 5000 

19.2714 

16.7714 

35.7561 

35.7561 

-10. 5373 

300 

2. 5000 

2. 5000 

19.2869 

16.7869 

35. 5510 

35. 5292 

-10.4417 

a 336. 35 

2. 5000 

2. 5000 

19. 5728 

17.0728 

31.9792 

30.7179 

-8.7772 

400 

2. 5000 

2. 5000 

20.0061 

17. 5061 

27.2883 

25.6189 

-6.6626 

500 

2. 5000 

2. 5000 

20. 5640 

18. 0640 

22.3306 

20. 2476 

-4.4500 

600 

2. 5000 

2. 5000 

21.0198 

18. 5198 

19.0255 

16.6804 

-2.9924 

700 

2. 5000 

2. 5000 

21.4052 

18.9052 

16.6647 

14. 1398 

-1.9629 

800 

2. 5000 

2. 5000 

21.7390 

19.2390 

14.8941 

12.2373 

-1. 1993 

900 

2. 5000 

2. 5000 

22.0334 

19. 5334 

13. 5170 

10.7568 

-.6120 

1000 

2. 5000 

2. 5000 

22.2968 

19.7968 

12.4153 

9. 5687 

-. 1474 

1030. 90 

2. 5000 

; 2. 5000 

22.3729 

19.8729 

12. 1181 

9.2471 

-.0231 

1100 

2. 5000 

! 2. 5000 

22. 5351 

20.0351 

11. 5139 

8. 5907 

.2281 

1200 

2. 5001 

: 2.5000 

22.7527 

20. 2526 

10.7628 

7.7677 

. 5370 

1300 

2. 5003 

2. 5000 

22.9528 

20.4528 

10. 1272 

7.0618 

.7946 

1400 

2. 5008 

2. 5001 

23. 1381 

20.6380 

9. 5824 

6.4457 

1.0119 

1500 

2. 5018 

2. 5001 

23.3107 

20.8105 

9. 1103 

5. 8994 

1. 1968 

1600 

2. 5034 

2. 5003 

23.4722 

20.9719 

8.6974 

5. 4082 

1.3552 

1700 

2. 5060 

2. 5005 

23.6240 

21.1235 

8.3331 

4. 9606 

1.4917 

1800 

2. 5099 

2. 5009 

23.7673 

21.2664 

8.0095 

4. 5478 

1.6097 

1900 

2. 5154 

2. 5016 

23.9032 

21.4016 

7.7201 

4. 1627 

1.7120 

2000 

2. 5228 

2. 5024 

24. 0324 

21. 5300 

7.4601 

3.7999 

1. 8007 

2100 

2. 5324 

2. 5036 

24. 1557 

21.6521 

7.2252 

3.4548 

1.8775 

2200 

2. 5446 

2. 5052 

24.2738 

21.7686 

7.0121 



2300 

2. 5596 

2. 5072 

24.3872 

21.8800 

6.8182 



2400 

2. 5775 

2. 5098 

24.4965 

21.9868 

6.6411 



2500 

2. 5988 

2. 5129 

24.6022 

22.0893 

6.4790 




a The reference element is crystal potassium below 336. 35° K and liquid potassium 
above 336. 35° K. 



TABLE m. - THERMODYNAMIC PROPERTIES OF IDEAL POTASSIUM GAS (Kg) 


Temper- 

ature, 

°K 

Heat cap- 
acity for 
standard 
state at 
constant 
pressure, 
C°/R 

Sensible 
enthalpy 
for stand- 
ard state 
at T° K, 

n T " H 0 
RT 

Entropy 

for 

standard 

state, 

S°/R 

Sensible 
free en- 
ergy for 
standard 
state at 
T° K, 

-K - H o) 

Sum of 
sensible 
enthalpy 
at T°K 
and 

chemical 
energy at 
0° K for 
standard 
state, 
H°/RT 

Formation from assigned 
reference state 

Enthalpy 
change 
at T° K, 
(AH°) f /RT 

Log of 
equilib- 
rium 
constant, 
log 

RT 

100 

4.3895 

3. 9898 

25. 1220 

21. 1321 

142.6636 

155. 5362 

-60.3574 

200 

4. 5137 

4.2285 

28.2136 

23.9851 

73. 5654 

76.8817 

-26.7246 

298. 15 

4. 5579 

4. 3303 

30.0248 

25.6945 

50. 8417 

50. 8417 

-15.7969 

300 

4. 5585 

4.3317 

30. 0530 

25.7213 

50 . 5563 

50. 5127 

-15.6610 

a 336. 35 

4. 5709 

4.3569 

30. 5751 

26.2182 

45. 5859 

43.0634 

-13.2975 

400 

4. 5905 

4. 3925 

31.3689 

26.9764 

39.0610 

35.7222 

-10.3406 

500 

4.6188 

4.4350 

32.3964 

27. 9614 

32. 1697 

28. 0037 

-7.2671 

600 

4.6455 

4. 4679 

33. 2408 

28.7730 

27. 5802 

22.8898 

-5.2585 

700 

4.6715 

4.4951 

33.9589 

29.4638 

24. 3056 

19.2558 

-3.8510 

800 

4.6972 

4. 5188 

34. 5844 

30.06 56 

21.8530 

16. 5393 

-2.8149 

900 

4.7226 

4. 5400 

35. 1391 

30. 5991 

19.9482 

14. 4277 

-2.0239 

1000 

4.7479 

4. 5595 

35.6380 

31.0785 

18.4269 

12.7337 

-1.4031 

1030. 90 

4.7557 

4. 5653 

35.7826 

31.2173 

18.0170 

12.2749 

-1.2379 

1100 

4.7730 

4. 5778 

36.0917 

31. 5139 

17. 1845 

11.3381 

-.9053 

1200 

4.7982 

4. 5951 

36. 5081 

31. 9130 

16. 1513 

10. 1612 

-.4994 

1300 

4. 8233 

4.6117 

36.8931 

32.2814 

15.2789 

9. 1481 

-. 1639 

1400 

4. 8483 

4.6277 

37. 2515 

32.6238 

14. 5330 

8.2594 

. 1161 

1500 

4. 8734 

4.6432 

37. 5868 

32.9436 

13.8882 

7.4664 

.3516 

1600 

4.8984 

4.6584 

37.9021 

33.2437 

13.3255 

6.7473 

.5507 

1700 

4. 9234 

4.6733 

38. 1999 

33. 5266 

12. 8305 

6.0856 

.7197 

1800 

4. 9484 

4.6878 

38.4820 

33.7941 

12. 3919 

5. 4686 

.8631 

1900 

4. 9733 

4.7022 

38.7502 

34. 0480 

12. 0008 

4. 8860 

.9846 

2000 

4. 9983 

4.7164 

39.0059 

34.2895 

11.6501 

4. 3297 

1.0873 

2100 

5. 0233 

4.7304 

39.2504 

34. 5200 

11.3339 

3.7931 

1. 1734 

2200 

5. 0482 

4.7443 

39.4847 

34.7404 

11.0476 



2300 

5. 0732 

4.7580 

39.7096 

34.9516 

10.7873 



2400 

5. 0982 

4.7717 

39.9260 

35. 1543 

10. 5498 



2500 

5. 1231 

4.7853 

40. 1347 

35.3494 

10.3322 




^he reference element is crystal potassium below 336. 35° K and liquid potassium 
above 336. 35° K. 




Liquid 


The following equations for the thermodynamic properties of the liquid in the range of 
338° to 1070° K were taken from reference 13: 


(h°) - (H°) = -231. 40909 + 8. 8849061 T - 2. 2831074xl0" 3 T 2 + 9. 792447X10' 7 T 3 (1) 
l c 


(C°) = 8. 8849061 - 4. 5662147x10 3 T + 2. 9377334x10“ 6 T 


\-6 rr .2 


( 2 ) 


- 20. 458252 log T - 4. 5662147X10" 3 T + 1. 4688667X10 -6 - 32. 310079 (3) 


The units for heat capacity and entropy are calories per mole per °K (4. 184 J/ 
(mole)(°K)). Equation (1) agrees well with measurements from 1040° to 1329° K 
(refs. 2 to 4). 


Monomer and Dimer (Ideal Gases) 

The methods used to compute the thermodynamic properties of the ideal monomer 

and the ideal dimer are the methods for monatomic and diatomic molecules described in 

reference 14. The energy levels for atomic potassium were taken from reference 15. 

The molecular constants for the dimer were taken from reference 16 and are shown 

in table IV. The chemical energy of the monomer in its standard state can 

' 'monomer 


TABLE IV. - MOLECULAR CONSTANTS 3 FOR DIATOMIC POTASSIUM 


Molecular weight, g/mole 

Symmetry number 

Electronic state 

Statistical weight 

Zero- order vibrational frequency, & e , cm -1 

Anharmonicity constant, co x , cm * 

e e 

Rotational constant, B , cm" 1 

e - 1 

Vibration- rotation interaction constant, a , cm 

e -1 

Spectroscopic constant for rotational stretching, D e , cm 
a Constants from Herzberg (ref. 16). 


78.204 

2 

*2 

1 

92.64 

0.354 

0.05622 
0. 000219 
8. 34X10' 8 
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be computed from the relation 


K) 


monomer 


= (H°) -(aH° 98 ) + a(h° 98 -H°) = 19 703.3 cal/mole (82 438.6 J/mole) 


where 


A ( H 298 " H o)„ ~ ( H 298 " H o) monomer " ( H 298 “ H o), 


The chemical energy of the crystal in its standard state (Hq) has been shown to equal 

-1693. 0 calories per mole (-7083. 5 J/mole). The average value of (AHggg) is 

21 184. 6 calories per mole (88 636. 4 J/mole) from table V. This value is the average of 
the apparent values resulting from the thermodynamic analysis presented herein. 


TABLE V. - APPARENT VALUES OF 
ENTHALPY CHANGE ON EQUILIBRIUM 
VAPORIZATION TO 1 MOLE OF REAL 
MONOMER AT 298. 15° K AS FUNC- 
TION OF TEMPERATURE 


Temperature, 

T, 

°K 

Apparent values 

a of(AH° 98 ) v 

cal/mole 

J/mole 

900 

21 256.2 

88 935.9 

950 

21 239.0 

88 864.0 

1000 

21 223.0 

88 797. 0 

1030.9 

21 213.6 

88 757.7 

1100 

21 194. 1 

88 676. 1 

1150 

21 181.0 

88 621.3 

1200 

21 168.6 

88 569.4 

1250 

21 156.7 

88 519.6 

1300 

21 144.9 

88 470. 3 

1350 

21132.9 

88 420. 1 

1400 

21 120. 5 

88 368.2 


a Average (AHggg) = 21 184. 6 cal/mole; 

= 88 636. 4 J/mole 


Spread =136 cal (569 J) 
= 0.65 percent 
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From reference 17, 


(ah£) = -D° = (h°) ~ 2 f H n) = -11 850 cal/mole; -49 580 J/mole 

' u/ dimer ' u/ dimer v 'monomer 

Therefore, 


HS) = -11 850 + 2(19 703.3) = 27 556.6 cal/mole; 115 296.8 J/mole 
k u ' dimer 


Heats of Formation and Equilibrium Constants 


Values for heats of formation of the ideal species from K (crystal and liquid) 

^i. e. , ^AH^ ^are given in tables I to in. For Kg (ideal gas) at 298. 15° K, for example, 

( iH 298. 15 ) (j dimer = ( H 298. l 5 ) dimer ' 2 ( H 298. 15> c 


The logarithms of the equilibrium constants for the type of reaction just discussed 
are also listed in tables I to m. The equilibrium constant for formation from the 
assigned reference element is obtained from the standard free- energy change ^AF^,j , 
for example, * 


(af°) = (f°) 

v 1/ f, dimer v i/ dimer 



cond 


by the following equation: 


( af t) 

log X (4) 

1 2. 3025851 RT 

Physical Constants 

The atomic weight, 39. 102 grams per mole, was taken from Cameron and Wichers 
(ref. 18). The values of the universal gas constant R, 1.98717 calories per mole per °K 
(8. 3143 J/(mole)(°K), the entropy constant, -3.664953, and the second radiation constant 
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hc/k, 1. 43879 were computed from the constants given in reference 19. 


SELECTION OF REAL-GAS DATA 

The initial data considered in the thermodynamic analysis are vapor pressures, sec- 
ond virial coefficients, heat of dissociation of dimer at 0° K, heats of vaporization, and 
volumes of vapor. Where it is appropriate, the manner of selecting these data is dis- 
cussed in the following paragraphs. 


Vapor Pressures 

References 2 to 4 present vapor pressures measured with great care from 1127° 
to 157 5° K. Several sets of data in the range from 372° to 2078° K are also available. 
Each set was combined with the data of references 2 to 4 and fitted by the least- squares 
technique to a three- constant equation. Four of the sets did not increase the standard 
deviation of the data of references 2 to 4 appreciably. Therefore, these four sets 
(refs. 9 to 11 and Bowles' data) were pooled with the data of references 2 to 4 and yielded 
an overall standard deviation of 1. 43 percent. 

The data of reference 20 were not used because of the uncertain corrections made to 
the raw data. The data of Buck and Pauly (ref. 21) and of Rigney et al. (ref. 22) are con- 
sidered acceptable but were not available when the analysis was completed. 


Temperature, T, °K 

2000 1000 666.7 500 400 333.3 



Figure 1. - Smoothed experimental vapor pressures above liquid potassium (336. 35° to 2100° K). 


9 



TABLE VI. - BOILING POINTS OF POTASSIUM FROM 
l(f 8 TO 10 2 UNITS OF PRESSURE REFERRED TO 
STANDARD PRESSURE 


Pressure, 

log(P/P 0 ) 

Boiling points, a 
°K 

Pressure, 

log(P/P 0 ) 

Boiling noints, a 
°K 

-8 

361.21 

-2 

700. 36 

-7 

392.64 

-1 

833. 18 

-6 

430. 17 

0 

1030. 90 

-5 

475.78 

1 

1359.35 

-4 

532.41 

2 

2027.78 

-3 

604.71 




a From eq. (5) 



-4 625.63 + 6 5981? _ 0 7006 43 log T 
T 


The equation fitted to the accepted vapor pressures referred to the standard pressure 
P Q of 1 atmosphere (1. 01325X10^ N/m 2 ) is 



4625 ‘ 3 + 6. 59817 - 0. 700643 log T 
T 


(5) 


Figure 1 and table VI present the smoothed experimental vapor pressures of liquid potas- 
sium from this equation as a function of the boiling points. 


Heat of Dissociation of Dimer at 0° K and Second Virial Coefficients 

The heat of dissociation of the dimer Dq at 0° K and the second virial coefficients 
can be related by the following expressions 

-B = RTcs^jirnerization ^ re ^' 


( iF °T) dlmerization ' ^ ln 23) 


and the identity 


(af°) = a(f°-h°) -d° 

'dimerization dimerization 


10 



where 


Therefore, 


Dr 


•( ah oL 


dimerization 



-B = RT exp 


- a(f° - H°) 

v 'dimerization 

RT 


( 6 ) 


Equation (6) was used to calculate values of the second virial coefficient B from 
Dq = 11 850 calories per mole (49 580 J/mole) (ref. 17) and values of Fq - Hq for the 
ideal monomer and dimer. These values of B may be compared with the two latest sets 
in the literature, those from Davies et al. (ref. 6) and those from Ewing et al. (ref. 2 
to 4). 

Davies et al. calculated values of B from two estimates of the potential energy 
curve for the interaction of ground- state K atoms, one from the Rydberg potential^ and 
the other from an adjusted Rydberg potential. Ewing et al. obtained values of B graph- 
ically from their PVT data. 

The following table shows that the values of B based on Dq (By, cm^/mole) seem 
to favor the average of Davies' two sets. 


Temper- 

Second virial coefficient -B, cc/mole, from - 


ature, 

T, 

Eq. (6) for a 
dissociation 
energy D° of 
11 850 cal/mole 


Ref. 6 


Ref. 2 

°K 

By Rydberg 
potential 

By adjusted 
Rydberg po- 
tential 

Average 


1100 

3900 

3550 

4390 

3970 

4824 

1200 

2620 

2320 

2890 

2605 

3276 

1300 

1860 

1630 

2040 

1835 

2376 

1400 

1420 

1250 

1570 

1410 

1815 

1500 

1111 

960 

1210 

1085 

1443 

1600 

906 

770 

990 

880 

1186 


*The Rydberg potential is "the most realistic three- parameter diatomic potential 
energy function that is available at present" (ref. 24). It uses spectroscopic values of 
Dq, equilibrium internuclear separations, and vibrational frequencies. 
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Therefore, a curve fit of the average values of B was used in the analysis. The equa- 
tion selected was of the following form: 

log(-B) = — + b + c log T (7) 

T 

The method of least squares yielded the constants 

a = 2943.72 
b = -2. 52548 
c = 1. 13330 

Heat of Vaporization of Liquid and Volume of Vapor 

Achener (ref. 9) reported direct measurements of the heat of vaporization of liquid 
potassium in the range 815° to 1315° K. Ewing et al. (refs. 2 to 4) have measured the 
volumes of the saturated vapor in the range 1127° to 1575° K. 

DEVELOPMENT OF VIRIAL EQUATION OF STATE 

If it is assumed that potassium vapor exists as an ideal monomer, it can be shown 
that third-law heats of vaporization at some reference temperature vary appreciably with 
the temperature at which the vapor pressure was measured. For example, substituting 
the fitted vapor pressures from equation (5) and the values for sensible free- energy 
changes for vaporization to the ideal monomer into the third-law expression 

( AH 29s) v - - RT >" p * a ( f T - H 29 s) v 

yields heats of vaporization of 20 850, 20 557, 20 010, and 19 146 calories per mole 
(87 236, 86 010, 83 722, and 80 107 J/mole) at 500°, 1000°, 1500°, and 2000° K, respec- 
tively. This nonideality can be explained by assuming either that the vapor consists of 
more than one molecular species or that it is a monatomic imperfect gas. 

Preliminary analysis in terms of an ideal mixture of monomer and dimer did not 
correlate the experimental vapor pressures satisfactorily. Since the properties of the 
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higher polymers of potassium are not known, the mixture approach was not pursued any 
further. Therefore, the vapor was treated as an imperfect gas. 


Formulation of Third and Fourth Virial Coefficients 

The virial equation of state, to the number of terms used in the present work, 

PV = RT + BP + CP 2 + DP 3 (8) 

describes the deviations from ideality in terms of the virial coefficients B, C, and D, 
which are functions of T. 

It is necessary to express B, C, and D as explicit functions of T in order to use 
the equation of state (eq. (8)) directly. Unfortunately, an explicit function based on a 
realistic potential (ref. 6) is available only for B (i. e. , eq. (7)). One way to overcome 
this difficulty is to assume relations of the following type: 


c4(-B) m 

d 2 (-B) q 


(9a) 


(9b) 


and to substitute these relations into equation (8). Then the constants d^, d 2 , m, n, q, 
and r could be fitted to the PVT data (refs. 2 to 4). 

The preceding empirical approach was improved by deriving values of m, n, q, and 
r from two less realistic potentials, the Lennard-Jones 12-6 (ref. 7) and the hard- 
sphere (ref. 8) potentials. 

In this way, only the scale factors dj and d 2 were fitted to the PVT data. The set 
of constants derived from the Lennard-Jones potential is m = 8. 5561, n = -10. 056, 
q = 12. 359, and r = -13. 093. For the hard-sphere potential, the constants are m = 2, 
n = 1, q = 3, and r = 2. For this potential, equations (9a) and (9b) simplify to 


C = 



(10a) 
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(10b) 


D = 



The hard- sphere relation reproduces the experimental volumes (refs. 2 to 4) some- 
what better than the Lennard- Jones relation (0. 75 percent standard deviation against 
2. 20 percent). Either potential provides a good representation of the experimental heats 
of vaporization (ref. 9) (3. 11 percent standard deviation against 3. 25 percent). 

Table VII compares the two models with respect to fugacity- pres sure ratio f/P, heat 
capacity Cp, and isentropic expansion index y g . Figure 2 shows the variation in Cp. 

From 1000° to 2000° K, either potential yields a gradual variation for f/P, C , and y . 

P 5 

However, between 800° and 1000° K, the Lennard-Jones potential shows a sharp decrease 

in C and a sharp increase in y . 

P s 

TABLE vn. - COMPARISON OF THERMODYNAMIC PROPERTIES FROM LENNARD-JONES 12-6 

AND HARD- SPHERE POTENTIALS 


Temper- 


Lennard-Jones potential 



Hard- sphere potential 

ature, 

T, 

Fugacity- 

Real- gas heat capacity at 

Isentropic 

Fugacity- 

Real- gas heat capacity at 

°K 

pressure 

constant pressure, C 

expansion 

pressure 

constant pressure, C 

ratio, 

f/P 



index, 

ratio, 

f/P 




cal 

J 

cal 

J 


(mole) (°K) 

(mole) (°K) 

^s’ 

/ a In p\ 
\a In P/ S 

(mole)(°K) 

(mole) (°K) 



336. 35 





0. 99995 

5.012 

20. 970 





500 

.9982 

5.627 

23. 543 





800 

0. 9307 

39.37 

164.72 

1.095 

.9744 

8.679 

36.313 

900 

.9564 

21.70 

90.79 

1.201 

.9590 

9.748 

40. 786 

1000 

.9418 

11. 36 

47. 53 

1.392 

.9408 

10.671 

44. 647 

1500 

.8420 

11. 59 

48.49 

1.451 

.8336 

12.311 

51. 509 

2000 

.7569 

10. 96 

45.86 

1.389 

.7423 

10.683 

44.698 


Isentropic 

expansion 

index, 

v 

/a In p\ 
\3 ln P/ s 

1.659 
1. 583 
1.448 
1.437 
1.434 
1.440 
1.371 


On the other hand, the hard- sphere relation shows a gradual variation in all three 
properties over the entire range of temperatures. Therefore, the hard-sphere relation 
probably gives a better description of the properties of potassium vapor. 

This conclusion does not imply that potassium is a hard-sphere gas. It may indicate 
that, when the Rydberg potential curves for higher interactions are eventually determined, 
the expressions for C and D will show dependencies on B which are similar to those 
for a hard-sphere gas. 
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Figure 2. - Comparison of calculated heat capacities of saturated vapor. 


Hard-Sphere Expressions Used for Higher Virial Coefficients 


Equation (7) derived from Davies' second virial coefficients can be expressed as 


By = -10 b 10 a / T T c cm 3 /mole 


( 11 ) 


Expressing B in calories per mole per atmosphere (4. 12929x10 5 J/(mole)(N/m 3 )) 
yields 


B = Cj 10 a / T T c 


(11a) 


where 


c x = -0. 024217256xl0 b 


Substituting equations (10) into equation (8) results in 
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( 12 ) 


p2p2 p3p3 

PV = RT + BP + k 1 B-f— + k 9 £_J_ 

1 rri 6 n 


Further substitution of equation (11a) into equation (12) yields 


PV = RT + Cj 10 a / T T C P + c 2 10 2a / T t 2c- 1 P 2 + c 3 10 3a / T t 3c “ 2 P 3 


(13) 


where 


c 2 = kjc 2 (14a) 

c 3 = Vl (14b) 

With the values of a, b, and c (and therefore Cj) determined from Davies' second virial 
coefficients, the constants c 2 and c 3 were fitted to the PVT data. The curve-fitting 
process yielded 


c 2 = -2. 55216x10“ 9 


and 


c 3 = 2. 41943x10 13 


Fitting equation (8) yields 


c x = -7.22185xl0“ 5 


Therefore, 


k x = -0. 489340 


and 


k 2 - -0.642345 
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Thus, the final recommended equation of state becomes 

PV = RT- 7. 22185x10" 5 10 2943 ' 72 / T p 1, 13330 p_ 2. 55216X10' 9 10 5887 - 44 / T p 1 * 26660 p 2 

+ 2. 41943X10" 13 10 8831 - 16 / T t 1- 39990 P 3 (15) 


Consistency Test for Equation of State 

It is necessary to calculate the heat of vaporization at 298. 15° K as a function of tem- 
perature in order to check the consistency of the calculated properties from equation (15). 
The basic equation for a real gas obeying equation (8) is 


(F - F°) t = RT In P + BP + 


(16) 


which is derived in appendix B. The relative free energy (F - F°),p is the constant- 
temperature free energy of the real gas relative to the free energy of the ideal gas. 
The quantity Q, defined in reference 5 as 


Q ■ -A(F° - H° 98 ) - RT In P 


(17) 


can be simply related to the virial coefficients in the following manner. The ideal free 
energy of vaporization is defined as 


(AF°) . (F°) 


monomer 


( f t) 


cond 


(18) 


On the saturation line, equilibrium demands 


(Fiji) = (Frp) 

monomer cond 


(19) 


If it is assumed that the condensed phase is in its standard state, 


(F T ) 


monomer 


= ( f t) 


cond 


( 20 ) 
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Substituting equation (20) into equation (18) yields 


( AF t) v = < F ° - F >T, monomer 
and substituting equation (21) into equation (17) yields 

« = < F - f °>T, monomer ' RT ln P + (aH° 98 ) 


( 21 ) 


( 22 ) 


Substituting equation (16) into equation (22) gives 


( aH 29a) =Q-BP- 


CP 2 DP 3 


(23) 


Using equations (10) results in 


k l B 2 P 2 k 2 B 3 P 3 


K 98 ) v = «-BP -i T 3 , 


(24) 


Since kj and kg have already been determined by fitting the PVT data (1127° to 
1575° K) to equation (12), it is a simple matter to compute apparent values of ^AHggg) . 

The choice of a reasonable interval for averaging ^AHggg) must be a compromise 

because there is only a small interval (1100° to 1300° K) that includes experimental 
values of both the vapor volume and the heat capacity of the liquid. A representative av- 
erage probably could be obtained by a small extrapolation of volumes to 900° K (i. e. , 
using eq. (13)) and a small extrapolation of the cubic equation for heat capacity (eq. (2)) to 
1400° K. The interval from 900° to 1400° K was thus chosen for averaging ^AHggg) . In 

this way, from 900° to 1100° K, experimental enthalpies were combined with extrapolated 
volumes, and, from 1300° to 1400° K, extrapolated enthalpies were combined with exper- 
imental volumes. 

Table V shows the variation in apparent values of (AHggg) as a function of temper- 
ature. The maximum difference is 0. 65 percent, with an average value of 21 184. 6 cal- 
ories per mole (88 636. 4 j/mole). 


REAL-GAS PROPERTIES ON SATURATION LINE 

Table VTO presents calculated real-gas properties on the saturation line at assigned 
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TABLE Vm. - THERMODYNAMIC PROPERTIES OF REAL MONOMER GAS ON SATURATION LINE 


Temper- 

Pressure 

Volume, 

Real- gas 

Enthalpy 

Enthalpy 

Enthalpy 

Entropy 

Entropy 

Entropy 

Isen- 

ature, 

referred to 

cm 3 /mole 

heat ca- 

(crystal 

change on 

of real 

for stan- 

change on 

of real 

tropic 

T, 

standard 


pacity at 

phase). 

equilibrium 

gas at 

dard 

equilibrium 

gas at 

expan- 

°K 

pressure. 


constant 

h£(c) 

vaporization 

T °K 

state 

vaporization 

T °K 

sion 


p/p 0 


pressure 

RT 

to 1 mole of 

(gas 

(crystal 

to 1 mole of 

(gas 

index, 




(gas 


real monomer 

phase), 

phase). 

real monomer 

phase), 

^s 




phase). 


at T °K, 

H T (g) 

s£(c)/R 

at T °K, 

S T (g) 





C p (g)/R 


(H t ) /RT 

V 

RT 


(as t ) /r 

V 

R 


a 336. 35 

1 . 18949x10 9 

2. 32015X10 13 

2. 5219 

1.2612 

31.0561 

32. 3174 

9.0651 

31.0561 

32.3174 

1.6592 

400. 00 

1. 62548x10 7 

2. 01864X10* 1 

2. 5850 

1.6694 

25. 9050 

27. 5744 

9.7284 

25. 9050 

27. 5744 

1.6403 

500. 00 

2. 856 92x10" 5 

1. 43345x10° 

2.8315 

2.0830 

20. 4411 

22. 5242 

10. 5629 

20.4411 

22. 5242 

1. 5832 

600. 00 

8. 75583x10 ^ 

5. 58893X10 7 

3.2585 

2.3452 

16. 7689 

19.1141 

11.2296 

16.7689 

19.1141 

1. 5198 

700. 00 

9.92504X10' 3 

5. 70606X10 6 

3. 7976 

2. 5249 

14.1122 

16.6371 

11.7851 

14.1122 

16.6371 

1.4741 

800. 00 

6.05489X10' 2 

1. 05570x10° 

4. 3676 

2.6569 

12. 0866 

14.7435 

12.2633 

12.0866 

14.7435 

1.4483 

900. 00 

2. 44746x10” 1 

2. 88902X10 5 

4. 9055 

2.7602 

10. 4818 

13.2420 

12.6858 

10. 4818 

13.2420 

1.4368 

1000. 00 

7.42360x10”* 

1.03634X10 5 

5.3699 

2.8466 

9.1731 

12.0196 

13.0675 

9.1731 

12.0196 

1.4341 

b 1030.90 

1.00001X10° 

7. 87392x10** 

5. 4945 

2.8710 

8.8155 

11.6865 

13. 1789 

8.8155 

11.6865 

1.4343 

1100.00 

1.82863X10° 

4. 51586X10 4 

5.7376 

2.9232 

8. 0817 

11.0050 

13. 4191 

8.0817 

11.0050 

1.4358 

1200.00 

3.85545x10° 

2.27353X10 4 

5. 9992 

2.9950 

7.1554 

10.1504 ! 

13.7483 

7.1554 

10.1504 

1.4390 

1300. 00 

7.21497x10° 

1.27796X10 4 

6.1561 

3.0654 

6. 3572 

9. 4226 

14.0612 

6.3572 

9.4226 

1.4416 

1400. 00 

1.22982X10 1 

7. 82923X10 3 

6.2171 

3.1368 

5.6605 

8.7973 

14.3624 ■ 

5.6605 

8.7973 

1.4421 

1500. 00 

1. 94587X10 1 

5. 13751X10 3 

6. 1954 

3.2109 

5. 0450 

8.2559 

14.6554 

5. 0450 

8.2559 

1.4396 

1600.00 

2. 89875X10 1 

3. 56476X10 3 

6.1065 

3.2891 

4. 4949 

7.7840 ' 

14.9433 

4. 4949 

7.7840 

1.4334 

1700. 00 

4. 10983X10 1 

2. 5902 lx 10 3 

5. 9662 

3. 3725 

3.9975 

7.3700 

15. 2286 

3.9975 

7. 3700 

1.4233 

1800. 00 

5. 59237X10 1 

1. 95624XI0 3 

5.7893 | 

3.4617 ; 

3. 5427 

7.0043 | 

15. 5130 | 

3. 5427 

7. 0043 

1. 4094 

1900. 00 

7. 35185X10 1 

1. 52666X10 3 

5. 5888 

3.5574 

3.1219 

6.6793 

15. 7985 

3. 1219 

6.6793 1 

1.3919 

2000. 00 

9. 38677x10* 

1.22 537X10 3 

5.3761 

3.6602 

2.7282 

6.3884 

16.0863 

2.7282 

6.3884 

1.3713 

2100.00 

1. 16897X10 2 

1. 00773X10 3 

5. 1603 

3. 7704 

2.3557 

6.1262 

16. 3777 

2.3557 

6.1262 

1.3482 


Temper- 

Fugacity- 

Second virial 

Third virial 

Fourth virial 

Compress- 

Density, 

Volume- 

Volume- 

ature, 

pressure 

coefficient, 

coefficient, 

coefficient, 

ibility 

g/cc 

temperature 

pressure 

T, 

ratio. 

B v , 

<V 

D v , 

factor, 


coefficient, 

coefficient, 

°K 

f/P 

cm 3 /mole 

3 2 

(cm /mole) 

3 V 3 

(cm /mole) 

Z 


(3V/3T) p 

cm^/(moleX°K) 

(av/a(p/p )) 

3 T 

cm /mole 

a 336. 35 

0. 999947 

-1. 230535x10° 

-5. 928597X10*° 

8.281390X10 30 

0. 999947 

1.685319x10' 12 

1.672289x10° 

-4.723930X10 20 

400. 00 

. 999700 

-6. 063725x10° 

-1.439605xl0 17 

9.909236X10 26 

. 999700 

1. 937050x10'*° 

1. 228322X10 7 

-3.008377X10* 6 

500. 00 

. 998166 

-2. 634548x10° 

-2.717544X10* 4 

8. 127188X10 22 

. 998162 

2.727819X10 -8 

7. 114748X10 4 

-1.217337X10* 2 

600. 00 

. 993986 

-3. 382215X10 5 

-4. 478853X10* 2 

1.719593X10 20 

. 993950 

6. 996327X10" 7 

2. 409950X10 3 

-1. 555273x10° 

700. 00 

.986148 

-8.020207X10 4 

-2. 518463x10** 

2. 292867x10*° 

.985962 

6. 852720x10**° 

2. 245802X10 2 

-1.412355X10 7 

800. 00 

.974375 

-2.781269X10 4 

-3.028659X10*° 

9. 562052X10* 6 

.973749 

3.703889X10' 5 

3. 929410x10* 

-4. 338637X10° 

900. 00 

. 959012 

-1.239835X10 4 

-6.018550x10° 

8.470583X10 15 

.957445 

1. 353467X10' 4 

1. 040741x10* 

-2. 989555X10 4 

1000.00 

. 940765 

-6. 578704X10 3 

-1.694513x10° 

1.26 5444X10 15 

.937573 

3.773103X10' 4 

3.662369x10° 

-3.614396X10 3 

b 1030. 90 

.934673 

-5. 557 559X10 3 

-1.209296x10° 

7.629116X10 14 

.930823 

4. 966012X10 -4 

2.770519x10° 

-2.054201X10 3 

1100.00 i 

. 920472 

-3. 957669X10 3 

-6. 132582X10 8 

2.755120X10 14 

.914877 

8.658820X10' 4 

1. 575692x10° 

-6. 560752X10 2 

1200. 00 

.898953 

-2.613701X10 3 

-2.674707X10 8 

7.935782X10* 3 

. 890196 

1.719878X10 -3 

7. 844476X10”* 

-1.611839X10 2 

1300. 00 

.876931 

-1.853325X10 3 

-1. 344830X10 8 

2.829285X10 13 

. 864364 

3.059732X10" 3 

4.352644x10”* 

-4. 988836x10* 

1400. 00 

. 854994 

-1. 388940X10 3 

-7. 553203X10 7 

1. 190892X10* 3 

.838149 

4. 99436 2xl0' 3 

2.622304x10”* 

-1.848296x10* 

1500. 00 

.833588 

-1.087586X10 3 

-4.631175X10 7 

5.717580X10 12 

.812206 

7.611078X10' 3 

1.683745X10”* 

-7. 895512x10° 

1600. 00 

. 813034 

-8. 822125X10 2 

-3.047268X10 7 

3.051699X10 12 ' 

.787068 

1. 096903X10' 2 

1. 136799X10”* 

-3.781573x10° 

1700.00 

.793541 

-7.365236X10 2 

-2. 123917X10 7 

1.775751X10 12 

.763134 

1. 509605x10” 2 

7.990948X10' 2 

-1.987930x10° 

1800. 00 

.775232 

-6. 296788X10 2 

-1. 552396x1 0 7 

1. 109633x10 12 , 

.740688 

1. 998837X10' 2 

5. 804938X10” 2 

-1. 128295x10° 

1900.00 

.758165 

-5.491063X10 2 

-1. 180530X10 7 

7.358534x10** 

.719906 

2 . 56 1276X10" 2 

4. 333468X10” 2 

-6.825667x10"* 

2000. 00 

. 742348 

-4. 86 8 9 56X10 2 

-9. 281880x10° 

5. 130144x10** 

.700879 

3. 19 1038x10” 2 

3.310014X10” 2 

-4.356520x10”* 

2100.00 

.727753 

-4.37886QX10 2 

-7. 507345x10° 

3.731686x10** 

.683629 

3.880191X10" 2 

2. 578182x10” 2 

-2.909786x10"* 


a Melting point. 
b Normal boiling point. 


& 
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temperatures from 336. 35° to 2100° K. These values are computed by the use of equa- 
tion (5) for the vapor pressures and of equation (15), the derived equation of state. De- 
tails of the calculations are given in the following paragraphs. 


Virial Coefficients 


The virial coefficients from equation (15) are in calories per mole per atmosphere 
/ -10 2 

squared (4. 07527x10 J/(mole)(N/m ) ) for C and in calories per mole per cubic at- 

' / _ 15 o /3 \ 

mosphere 14.02199x10 0 J/(mole)(N/m ) ) for D. The values given in table VIII have 

V 3 2 7 3 3 

been converted to units of (cm /mole) and (cm /mole) , respectively, in order to con- 
form to the usual virial equation: 


PV = j , B V , C V , D V 
RT V 


V 2 


The conversion is accomplished by the following equations from reference 25: 


B = B, 


C = 


- B 


RT 


D = 


Dy “ 3ByCy + 2By 
(RT) 2 


Volume, Entropy, Enthalpy, and Heat Capacity 

The volumes are calculated at assigned temperatures from equation (15) by use of 
the smoothed experimental vapor pressures. The fugacity- pressure ratio f/P is calcu- 
lated from equation (B9). 

The relative functions (S - S°) T , (H - H°),p, and (C p - C p j are calculated from 

equations (Bll), (B12), and (B13). Adding the relative functions to the corresponding 
functions of the ideal monomer yields the following functions for the real gas on the sat- 
uration line: 
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S T = (s - s°) T + s° 

(24a) 

H t = (H - H°) t + H° 

(24b) 

C n = ( C n “ C n) + C n 

P V P P/ T P 

(24c) 


It is necessary to calculate H^, for the ideal monomer as follows, in order to com- 
pute the function Hrp for the real monomer: 


( iH 29 8 )„ - ( H 298)„ - ( H 298) 


monomer 


' "( H 298) 


monomer 


Since - Hggg) has already been calculated for the ideal monomer (table II, p. 4), 


( h t 


t) 


monomer 


= ( H T - H ! 98 ) 


monomer 


( AH 298) t 


Enthalpy of Vaporization and Entropy of Vaporization 

The enthalpy of vaporization in calories per mole (4. 184 J/mole) is simply 

(AH t ) = (Hrp) - (h°) 

v monomer 'I 

The entropy of vaporization, in calories per mole per °K (4. 184 J/(mole)(°K)) is calcu- 
lated from the expression 


(A Hrp) 

V 

v T 

The use of this expression in the real-gas table implies the use of the apparent values of 
• However, this use is necessary to meet the requirement of a zero Gibbs free 

energy change (AFm) on the saturation line. 

1 v 


(as t ) 
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Density 


Since the volume V of the real gas is in cubic centimeters per mole of monomer 
(39. 102 g), the density is p = 39. 102/V in grams per cubic centimeter. 


Volume-Temperature and Volume-Pressure Coefficients 


The volume-temperature coefficient (9V/9T) p and the volume- pressure coefficient 
(dV/dP)rj, are needed to compute the isentropic expansion index (eq. (29)) and also the co- 
efficients of thermal expansion and compressibility as follows (ref. 23): 


a = 




P 


and 


K = - 




Solving equation (8) for 


V and differentiating with respect to P at constant T yield 


V = — + B+CP + DP 2 
P 

and 

(—) = + C + 2DP 

\0P/ T p 2 

Differentiating V with respect to T at constant P results in 

/dv\ = R + dB + p dC + p 2dD 
\0T/p P dT dT dT 


(25) 


( 26 ) 


Differentiating equations (10) gives 
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(27a) 


dC _ 2C dB _ C 
dT B dT T 

dD _ 3D dB 2D 
dT B dT T 

Substituting equations (27) into equation (26) yields 

/V\A _ R + dB + 2PC dB + 3P 2 D dB _ PC _ 2P 2 D 
\3T/p P dT B dT B dT T T 

B + d ln B (B + 2CP + 3DP 2 ) - - (C + 2 PD) 
\9T/p P dT T 


(27b) 


(28) 


Isentropic Expansion Index 


The isentropic expansion index is defined (ref. 1) as 


7 


s 


v Zap 
p W 


s 


(29) 


where 



Simple substitution of T, C , equation (25), and equation (28) into equation (29) 
yields the isentropic expansion index. 


COMPARISON OF CALCULATED REAL -GAS PROPERTIES WITH EXPERIMENTAL 

DATA AND ANOTHER ANALYSIS 

The real- gas properties presented compare favorably with the experimental data and 
the calculated heat capacity appears to show a more reasonable trend with temperature 
than the result of an empirical analysis (ref. 2) as shown in the following paragraphs. 
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Enthalpy of Vaporization and Entropy of Vaporization 


Figure 3 compares Achener’s (ref. 9) experimental enthalpies of vaporization with 
those computed herein and those computed by Ewing (ref. 2). Figure 4 is a similar com- 
parison for the entropies of vaporization. 



100 


o 


90 


70 — 
60 — 


50 



Figure 4. - Comparison of calculated and experimental entropies of vaporization. 
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The two sets of calculated heats of vaporization agree very well and differ from ex- 
periment by about 3. 5 percent (maximum). The entropies show a similar trend. 


Volume 

When the choice of a potential function for the higher virial coefficients was dis- 
cussed in the section DEVELOPMENT OF VIRIAL EQUATION OF STATE, the author ex- 
plained that the constants kj and k 2 were fitted to the volumes, vapor pressures, and 
virial coefficients. The heats of vaporization were not fitted because their formulation 
not only involves the virial coefficients but also their derivatives with respect to temper- 
ature (cf eq. (B12)). It is therefore not surprising that the volumes deviated less from 
experiment than did the heats of vaporization (0.75 percent against 3. 11 percent). 


Isent roplc Expansion Index 


Goldman (ref. 1) reported experimental expansion indexes y g varying from 1.33 

K, with an average of 1. 38. From the equation 
of state of the present work, y e 


to 1. 42 in the range from 943° to 1051° 


s 



Figure 5. - Comparison of calculated and experimental compressibility factors. 


values of 1. 434 at 1000 u K and 
1.436 at 1100° Kwere computed. 


Compressibility Factor 

Figure 5 compares Ewing's 
nine experimental compressibility 
factors Z (ref. 2) with those com- 
puted by Ewing et al. (ref. 2) and 
those computed herein. Some of 
the experimental compressibility 
factors are closer to Ewing's curve, 
some are closer to the present 
curve, and a few agree well with 
either curve. In view of the scatter 
of the experimental points, neither 
set of calculated compressibility 
factors appears to be superior. 
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Heat Capacity 


Figure 2 and table VII illustrate the agreement of the hard-sphere and Lennard- 

Jones models in predicting a gradual increase in C p of the real gas from 1100° to 

1500° K, followed by a gentle decline. On the other hand, Ewing’s analysis predicts a 

gradual rise to 1200° K, a gradual decline to 1500° K, and then a sharp increase. The 

agreement of both potential models in this region tends to confirm a gradual change in 

C in this temperature range. 

P 


CONCLUDING REMARKS 


The method of analysis used herein, fitting of pressure-volume-temperature (PVT) 
data, may be compared with an alternative method of Hicks in which volumes are traded 
off with free energies (PFT data) . 

The advantages of the PVT fit are that it preserves the experimental vapor pressures 
and the properties of the liquid (specific heat, enthalpy, entropy, and Gibbs free energy), 
and only slightly perturbs the experimental volumes. The disadvantage of this approach 
is that it does not guarantee a temperature-independent value of enthalpy change on equi- 
librium vaporization to 1 mole of real monomer at 298° K (AH^ggj (eq. (23)). However, 

this variation in (AHggg) is small compared with the variation obtained by assuming 

ideal behavior. For example, between 1000° and 1400° K, ^AHggg^ for the ideal gas 


varies from 20 557 to 20 119 calories per mole (86 010 to 84 178 J/mole), whereas for 
the real gas, it varies from 21 223 to 21 120 calories per mole (88 797 to 88 366 J/mole). 
The variation has thus been reduced by 75 percent. 

The method of Hicks for fitting pressure, free energy, and temperature does yield a 
constant ^AH^gg) , but applying it to the present set of initial data caused a 9. 5-percent 


deviation in volume and approximately a 15-percent deviation in vapor pressures below 
800° K. On the whole, therefore, the PVT fit does a better job than the PFT fit. 

The variation of apparent (AHggg) obtained through this analysis at temperatures 

from 1127° to 1575° K may be caused by 

1. Truncation of the virial expansion at the fourth coefficient 

2. Assumption that the higher coefficients are simple functions of the second virial 

coefficient. 

3. Discrepancies in the experimental data. 

The third and fourth virial coefficients of the present analysis depart markedly from 
hard- sphere behavior (ref. 7) in two ways: First, the scale factors differ as follows: 
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Fit to PVT data (present work): 


C y = 0. 0276 B^. 

D y = -4.454 By 

Theoretical hard-sphere model: 

C y = 0.625 By 
Dy = 0. 2869 By 

where By, Cy, and Dy are the second, third, and fourth virial coefficients, respec- 
tively. Second, and of more importance, in the hard- sphere model the virial coefficients 
are independent of temperature, whereas Cy and Dy of the present analysis depend 
strongly on temperature. This dependence is caused by the fact that B is a function of 
temperature (eq. (7)). 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, June 16, 1967, 

129-01-02-01-22. 
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APPENDIX A 


SYMBOLS 


constant 

D V 

fourth virial coefficient, 

* c% 

second virial coefficient, 

O X 

(cm /mole) 

cal/(mole)(atm); 

d 

constant 

J/(mole)(N/m^) 

F 

Gibbs free energy 

rotational constant, cm"^ 

f t 

Gibbs free energy of real gas 

second virial coefficient, 


at T° K, cal/mole; J/mole 

cm^/mole 

( f t - H S) 

sensible free energy for 

constant 

standard state at T° K, 

third virial coefficient, 


cal/mole; J/mole 

2 

cal/ (mole) (atm) ; 

F o 

Gibbs free energy at 0° K, 

9 ^ 

J/ (mole)(N/m‘ i ) 

cal/mole; J/mole 

real-gas heat capacity at 

AF° 

standard free- energy change 

constant pressure, 

(AFJ 

Gibbs free- energy change on 

cal/ (mole) (°K); 

V 

vaporization to 1 mole of 

J/(mole)(°K) 


real monomer at T° K, 

standard- state heat capacity 


cal/mole; J/mole 

at constant pressure, 

f 

9 

fugacity, atm; N/m 

cal/(mole)(°K); 

g 

gas phase 

J/(mole)(°K) 

H 

enthalpy 

third virial coefficient, 



3 2 
(cm /mole) 

h t 

enthalpy of real gas to T° K, 
cal/mole; J/mole 

constant 

h° 

sum of standard- state sensi- 

fourth virial coefficient, 

3 

1 

ble enthalpy at T° K and 

cal/(mole)(atm) ; 

9 3 


chemical energy at 0° K 

J/(mole)(N/ni ) 


for standard state, 

spectroscopic constant for 


cal/mole; J/mole 

rotational stretching, 

rr° ttO 

h t h o 

sensible enthalpy for stan- 

- 1 

cm 

dard state at T° K, 

dissociation energy at 0° K, 


cal/mole; J/mole 

cal/mole; J/mole 





- RT In P, 


AH° 



(ah t ) 

v 



hc/k 

$ 


3 %^ 


k 

l 

m 

n 

P 


chemical energy at 0° K for 
standard state, cal/mole; 
J/mole 

enthalpy change for formation 
of substance from element in 
atomic gas state at T° K, 
cal/mole; J/mole 

enthalpy change for formation 
of substance from assigned 
reference element at T° K, 
cal/mole; J/mole 

enthalpy change on equilibrium 
vaporization to 1 mole of 
real monomer at T° K, 
cal/mole; J/mole 

enthalpy change on equilibrium 
vaporization to 1 mole of 
ideal monomer at T° K, 
cal/mole; J/mole 

second radiation constant 

equilibrium constant for reac- 
tion of formation from ele- 
ment in atomic gas state 

equilibrium constant for reac- 
tion of formation from as- 
signed reference element 

Boltzmann constant 
liquid phase 
constant 
constant 

2 

absolute pressure, atm; N/m 

standard pressure, 1 atm; 
1.01325X10 5 N/m 2 


Q 


q 

R 


r 

S 


Sip 




T 

V 

Z 

O! 


(X . 


Yc 


“ A ( F t “ H 298) v 
cal/mole; J/mole 

constant 

universal gas constant, 

1.98717 cal/(mole)(°K), 
82.05601 (cm S )(atm)/(mole)(°K); 
8. 3143 J/(mole)(°K) 

constant 

entropy of real gas, 

cal/(mole)(°K); J/(mole)(°K) 

entropy of real gas at T° K, 
cal/(mole)(°K); J/(mole)(°K) 

entropy for standard state, 
cal/(mole)(°K); J/(mole)(°K) 

entropy change on equilibrium 
vaporization to 1 mole of real 
monomer at T° K, 
cal/(mole)(°K); J/(mole)(°K) 

entropy change on equilibrium 
vaporization to 1 mole of 
ideal monomer at T° K, 
cal/mole; J/mole 

absolute temperature, °K 
3 

volume, cm /mole 

compressibility factor; PV/RT 

coefficient of thermal expansion, 
(°K) 

vibration-rotation interaction 
constant for diatomic mole- 
cule, cm - * 

isentropic expansion index, 

(9 In P/9 In p) g 
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K 

coefficient of compressi- 

1 2 '* 
bility, atm ; (N/ni ) 

f 

formation from assigned ref- 
erence element 

P 

o 

density of vapor, g/cm 

l 

liquid 

w e 

zero- order vibrational fre- 

P 

constant pressure change 


quency for diatomic mole- 

s 

constant entropy change 


cule, cm”* 

sat 

quantity at saturation 

(j) X <X> V 

ee’e y e 

anharmonicity constants for 
diatomic molecules, 

T 

constant temperature change 

Subscripts: 

c 

cond 

cm”* 

crystal 

condensed- phase property 

V 

0 

Superscript: 

o 

change on equilibrium vapor- 
ization to monomer 

quantity at 0° K 
standard state 
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APPENDIX B 


DERIVATION OF RELATIVE FUNCTIONS AND FUGACITY 
Relative Free Energy, (F - F°Fp 

For a system of constant composition, 

dF = V dP - S dT 


and 


I = V (Bl) 

T 

For 1 mole of ideal gas, V = RT/P, so that integration of equation (Bl) gives 

F = RT In P + F° (B2) 

where the integration constant F° depends only on the temperature and nature of the gas. 
For a real gas, the fugacity f is defined by the equation 

F = RT In f + F° (B3) 

As the actual pressure P approaches zero (i. e. , as the real gas approximates ideal be- 

havior), f/P approaches 1. 

To determine f, differentiate equation (B3) and utilize equation (Bl) to obtain 

) (B4) 

T 



= V = RT 


9 In f 
9P 



At a definite temperature, 


RT d In f = V dP (B5) 

where V is the actual molar volume. For the ideal gas 
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RT d In P = dP 


Subtracting equation (B6) from equation (B5) yields 


RT d In/— -) = (V - — |dP 


Integrating between the pressures zero and P gives 

P=P # ’ P 
P P=0 


RT 


In — 


/ 


PV - RT 


dP 


Because f/P approaches 1 at low pressures 


/ 


RT In f = RT In P + / — ~ RT dP 


(B6) 


(B7) 


(B8) 


The integral must be evaluated either graphically or analytically to evaluate the fu- 
gacity. Solving the equation of state for (PV - RT)/P and substituting this expression 
into equation (B8) give 


RT In f = RT In P 



(B + CP + DP 2 )dP 


or 


RT In f = RT In P + BP + 


CP 2 

2 



From equation (B3) 


RT In f = F - F° 

Substituting equation (BIO) into equation (B9) yields 


(B9) 


(BIO) 
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( 16 ) 


(F - F°) t = RT In P + BP + 

1 2 3 


which is equation (16) of the main text. 

Relative Entropy (S - S°)j 

For a system of constant composition, 


(— 1 = -S 

\oT 


Differentiating equation (16) yields 


S - S u = - 


9(F - F°) 


3T 


= -R In P 


JP 


dB p dC P 2 dD P 3 
dT dT 2 dT 3 


Relative Enthalpy (H - H°h 


Since F = H - TS, 


H - H° = (F - F°) + T(S - S°) 


= RT In P + BP + + 5*^ - t(r In P + — P + — ^ + ® ^ 

2 3 \ dT dT 2 dT 3 , 


= [B - T —V + (c - T — + (d - T — 


dT 


dT/ 2 


dT/ 3 


(BH) 


(B12) 
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Relative Heat Capacity (Cp - 

The relative heat capacity is simply obtained by differentiation of equation (B12) 


C - C° - 
C P C P 


0(H - H°) 

= -T 

/d 2 B\ p t / d 2 c\ P 2 ( /d 2 D\ P 3 

0T 

— X 

p 

VdT 2 / \dT 2 / 2 VdT 2 / 3 _ 


(B13) 
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